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ABSTRACT 


Hyperbolic paraboloid shells have good adoptability towards 
foiindations. The ease in construction and the advantage of the 
shape for confinement of the soil enable the use of the hyper- 
shells for foundations in lose soils. Unlike the roof structures, 
the foundation shells are subjected to very high intensities of 
loads with more critical secondary conditions. Therefore, the 
foundation shells are relatively thick and call for bending analy- 
sis, ''The primary object of the thesis is to develop a qualitative 
study of shallow, and thick hypar shells and then suggest a range 
of moment coefficients for veiy simple boundary conditions. The 
total strain energy of the shell, including the shear deformation 
and normal stress is considered and the governing equations are 
developed by variational principle. Deflection coefficients are . 
generated by Baleigh-Rits’s method for three different boundary 
conditions. The comparison of maximum deflection and bending 
moment of a simply supported shallow hypar shell Is made with thin 
and thick plate theory results. As the rise of the ^ell decreases, 
the deflection and bending moments appear to tend towards the those 
of simple plate. 

^4he deformation and bending behaviorir of the shell with 
respect to different parameters such as span, height and thickness 
ratios, (the thickness ratio is defined as the ratio of the thiclo- 
nese to short span, similarly the other ratios), is studied j!^\_The 



deflection and "bending moment coefficients are found to be very 
sensitive and change rapidly for thickness ratio varying from 
0.025 to 0.075. The variation of the coefficients damps and 
becomes asymptotic for the thickness ratios higher than 0.075. 

series of deflection and bending moment coefficients are pre- 

■■y ■ 

sented through a set of graphsv'^^ These coefficients can be used ; 
for design purpose of some idealised simple boundary conditions^ 
Similar coefficients for different type of boundary conditions ; 
can easily be generated, A design example given at the end of 
the text illustrates the use of the graphs and suggest a method- 
of design! ■ 



GHAPIER I 


1,1 IISTROIIUCTIOI: 

Theoretically, a shell may have any curved shape which is 
pleasing to the human feelings. However, the builders of the 
shells have been confined to a limited forms because of practical 
considerations. This last but very important factor speak heavily 
for the use of some ’anticlastic surfaces' which can be generated 
by straight lines and at the same time being clearly defined geo- 
metrically, j 1 j , Conoids, hyperboloids, and hype rbolicparabolo ids 
belong to this classification of shells. The last variety of 
shells are the only warped siirface whose equations are simple 
enough to permit stress cal culat ions, | 1 j , with comparative ease. 

Being a warped surface and easy to analyse, the hypor* 
shells have attracted many researchers to advance simplified ana- 
lytic procedure j2, 3l, Since most of the shells used, till very 
recently, have been for roof structures (i.e, where loading is 
nominal), therefore they have been assumed as thin shells'*'. 

This assumption was, in part responsible for the inhibited 
use of hyper shells as substructure, despite its clear superiority 
as load carrying member to 'plate’ foundation (i.e, radCt or mat 


* Hyperbolic paraboloids are commonly termed as hypar shells. 

+ A shell is said to he ’thin' if the thickness to minimum radius 
of curvature ratio can be neglected in comparison to tmity |4-l# 
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foundation). The nonavailability of a good theory for ’mode- 
rately thick’ shells has restricted the use of hypar shells for 
foundation structures to a very limited extent {5|. 

That the hypar shell foundation is far superior to 
ordinary mat foundation can be shown easily. The hyper shell 
foimdation (a typical of which is shovsm in Pig. 1,1) has all 
the advantages as that of a mat namely, 

(i) it provides an increased foundation area to develop 
minimtim contact pressure and maximum safety factor ggainst 
soil failure, 

(ii) settlement will be minimum if the compressible soil 
strata is at a shallow depth. Moreover, this will be pronounced 
to a greater degree in case of a hyper shell due to the developed 
skin friction (ELg. 1.2). This skin friction being inclined has 

a vertical component and thus offers resistance to settlement. 

In case of a mat it is almost Insignificant, 

(iii) If there are small localised weak spots, the foundation 
can provide adequate safety against that. 

(iv) Large hydrostatic lift can be resisted due to the heavy 
downward load. Furthermore, the continuity of the structure 
makes it possible to construct it soil flow tight. 

One of the draitoack with mat is that, the contact pressure 
distribution is nonuniform (Pig. 1.3) 16|. This cause large shear 
and bending stresses at the outer edges, particularly in case of 
a homogeneous, saturated clay. This defect can be obviated to a 
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large extent in a shell foundation, due to the mutual readjust- 
ment of skin friction forces. Again this more uniform pressure 
distribution preTents large concentration of bending and shear 
stresses at the outer edges and thus reducing the chances of 
’flipping up’ of the edges. A uniform pressure distribution along 
with skin friction makes its existence more welcome by lessening 
the chances of pxmching shear. 

With all these in perspective a necessity was felt for pro- 
viding a reasonably accurate procedure for design of such elements. 
And this work is only a beginning towards this goal without claim- 
ing to achieve a full proof design procedure. 

1.2 OBJEGO} MB SCOPE; 

This study is concerned with the bending analysis of shallow 
hypar shells with rectangular plan«-form and bounded by the charac- 
teristics. 

Since the motive force behind this work was to make available 
a set of data useful for actual design, many a simplifying assump- 
tions have been made, which might give a feeling of nausea to a 
theoretician, but it was made only when it was felt that no great 
sacrifice of acctiracy is involved# The price paid thus was well i 
worth considering the conse(iuent ease of design. 

following this guideline an approximate solution hy Eaheigh- 
mtE method was obtained. The results, for simply supported case, 
are presented in the form of graphs. The deflection function was 
ohosen to satisfy the boundary conditions of the problems. Shells 
with more difficult boundary conditions can be solved by Galerkin 



7 


method wherein only the geometric houndary conditions are satisfied, 
but it was not attempted here. 

Ihree boundary conditions were solved for, namely, 

(i) when all edges are simply supported, 

(ii) when two edges are clamped and two edges are simply 

suppo rted, 

(iii) when all edges are fixed, 

The first case was investigated in detail, and the design 

charts given are for this case only. The two other cases were 

solved to checlc the feasibility of the method. 

The simply supported case can be approximated in practice 
if, 

(i) there are construction joints along eg;,. and hf (Fig. 1.1), 

(ii) it is assumed that the edges ab, be, cd, da (Fig. 1.1) are 
supported by deep but narrow beams, 

(iii) it is assumed that the lines joining the columns forms a 
sort of simple support (which is in fact the assumption for flat 
slab design). 

1.3 SHELL iS oUHETIONS: 

The shells considered is shown in Fig. 1.4. They are doubly 
curved, shallow translational shells bounded by its characteristics 
and rectangular in planform. The middle surface of the shell, 
expressed in cartesian coordinates is given by: 

Z(x, y) = kxy (1.3.1) 

where k = c/ab 

and X and y are the characteristics of the hyperbola. 
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The shallowness is defined in the sense of Ylasovj7j. 

In this analysis besides the assmptions of the classical theory 
the following additional assumptions are made: 

(i) Asstimptions conceming G-eometry: The square of the partial 

derivatives of the surface function (1.5.1), Z* and Z and the 

X y 

products of these derivatives are negligible in comparison with 
unity. 

This assumption implies that the distance between two 
points on the surface can be replaced by its projected length. 

(ii) Assumptions concerning Internal Forces: Since the shell 
considered is shallow, the effects of "bending moments are well 
pronounced. 

(iii) Assumptions concerning Deformations: The effects of the 
tangential displacements u, v on the changes of curvature and 
twist of the surface are negligible. 

1.4 ”1111011 FOI^roiATIOlT: 

1.4.1 Introductory Differential Geometry: 

The surface (Fig. 1.4) is defined by 

Z = Z(x, y) (1.4.1) 

where, x, y, and z form a rectangular cartesian system. The 
radius vector to the middle surface is given by 

7 = xi + yl + (1.4.2) 


* For typing facility 8z/9x, az/8y etc. are written as z^, z^ 
respectively. A double suffix will denote a second order 
partial derivative. 
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r = xi + yj +• 2 (x, y) k 


(1.4.3) 


where i, j and k are the unit vectors along x, y, and z respectively 
The equation to the ..tangent vectors s and t to the surface can he 
deduced as: 


— ^ 6k •t- 
7+12 t: 


(1.4.4) 


The vectors sdx = ds and tdy = dt describe an elemental 

area on the shell surface. The element is not rectangular but 

its prog-ection on the xy plane is. The vectors s and t are not 

unit vectors. The unit vectors s, t and ^ (normal) direction 

are defined by e ,6+ and e_ , where, 

s' t ' 


s _ 


T + (6z/6x) k 

V 1+(8z!/0x)^ 


X J + (6z/6y) k 

« "«■ ", - * 

V 1+(8z/5y)^ 


(1.4.5) 


e = e„ X e . 
«; s t 


k -(6z/6x) i -(92/6y) 3 
■f 1 +( 0 z/ax)^ iri+(9z/0y)' 


and the angle os between s and t can be obtained as 


Cos c.> *= 


^ . t 


Is| |t 


(6z/6x) (Qz/6y) 

Yl+(dz/5x)^ f1+(az/6y)^ 


( 1 . 4 . 6 ) 
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and, 

sinw = ( . e . z/ 8 y) I, (1.4.7) 

iri+(6z/8x)^ '^1+(9z/9y)^ 

The element of surface is obtained as 

dA = If 1-r (c' 25 /cx)^ + (ez/Sy)^ dx dy (1.4.8) 

1.4.2 Pseudo-Stress Resultants; 

A representative small element of the shell surface is 
formed by two radial planes Tfshose horizontal lines are parallel 
to x-axis and the other radial planes of which the horizontal 
lines are parallel to y-axis, (Pigs. 1.6 and 1.7). 

The direct forces and T„ 4 . are taken as positive when 

xs yt 

it creates tension. Similarly, the inplane shears and Ty^ 
are positive when it creates tension in the diagonal direction of 
increasing values of x and y. and Ty^ are the transverse 

shears* 

The moments are positive when they produce tension in the 
top fibre. The positive twisting moments p 3 X)duce in plane shear 
stresses i such that the> result in diagonal tension along the in- 
creasing X and y direction. 

"the forces and moments are expressed per unit length. 

The force system in the shell element is reduced in terms 
of the force system in the plane element by resolving the vectors 

properly. 

The relation between forces and stresses are (Pig. 1.5)5 
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T 


sx 


‘sy 


T 






ty 


T 


tx 


= 


4'll/2 

-h/2 

+V2 

-'■ tgt »Kt +?)/Rt> 

-h/2 

+h/2 

I '^sa; 

-h/2 

+h/2 

■^_^/2 ^tt «Rs +’^)/V 

+h/ 2 

-V2 

+h/2 

+^)/Rg) de; 
-li/2 s B 


Similarly for moments we have 


and, 


+h/2 

+h/2 




+11/ 2 

-/ , 
-h/2 

+h/2 

/ 

-h/2 




de; 


ds: 


(1.4.9) 


(1.4.10) 


(1.4.11) 


( 1 . 4 . 12 ) 


In general T^y and are not equal though T and are. 
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This is due to the inequality of the radii of curvature in s 
and t directions. 


let 1 *^ and he the resultant forces and and 

are the net moments (per unit length) acting on the two :j&aces 
ahcd and hcef of the shell element (Pig. 1.6 and 1.7). Then, 


K = 
8 


T e+T e. + T_e^ 
sx s sy t s & 


or, 


I + (Qz/dx) _k j ^ I + (Sz/eyl_k ^ 


fl+CSz/Sx)^ 


“SX 


Vl+(8z/dy)^ 


sy 


+ ^ - ( 025/8x) i - (Qz/SiLl g. (1.4.15) 

lf1+(dz/5x)^ f^+{Qz/QJ)^ 


Again, writing down the resultant force per unit length 
in terms of the forces in the plane element we have. 


F dt 

s 


(T^^ i I + T ^2 k) dy 


XX 


xy 


or, 


T + T T + T k) 

'‘-^xx xy xz - — — 7772 

fl+(Bz/8y)^ 


(1.4.14) 


comparing (1.4.13) and (1.4.14) we have, 

rp ^ Vl+(8z/8y)j rp _ 

fi+(8z/ax)^ Y'i+(az/dx)'^ 


T = T 
xy sy 


( 8z/8y) 
fl+(0z/ey)^ 



(1.4.15) 


T. 


f 8 z/ 6 x) 


xz 


fl+(62/0x)' 


T + 
sx 


dz 

dy 


T + 
sy 


fl+( 02 / 6 x)‘ 


T 


s? 
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Writing "fl+Zy == K^. , 1+(5z/dx)^ = and 

Tl+CSz/dx)^ + (Sz/6y)^ = Q, the equations (1,4.15) can he 
rewritten as 

^xx ~ (^/^) " (VK^) (Sz/6x) 

^xy = I’gy - (1/V (1.4.16) 

^xz ’ ^VV ^^2/5x) Tg^ + (Bz/dy) + (l/K^) 

Solution of T , S and T gives 

o 2 »- ® ^ 

I = (K T - ( 8 z/ 8 x). I + I 9 |/ 8 £i j ) (e /Q) 

•^SX '^^y -^xx iy xy xz' ± 

I-, = (1/Q) (-(8z/ex)(ez/ay) ^ \y * 

(1.4.17) 

Similarly, 

Tyy = (VV’^ty - ((8®/9y)/V 

Tyjj = T^X ■ ((®z/®^)/V (1.4.18) 

Ty^ = (8z/5y) I^y + (8V8^) Ttl ■" 


and, 
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^ty ^ ^yy " ((Qz/9x)(ez/6y)/K^) 5?^^ 

+ ((82/6y)/K^) ) (Ky/Q) ^ 

®tx ^ (-C 625 /ex) ( 6 z/ey) Tyy + k:^ + (dz/dx)Ty^) (1.4.19] 

= -(Ky/Q) (( 6 z/ey) H- ( 6 z/ 0 x) - Ty^) | 

With the same argiiment eis before the moment relations are obtained 
as (from Pig. 1.7), 



In terms of the moments in the plane element, can be expressed 


H = 
s 


V"l+(6z/ey)‘ 




Comparing (1.4.20) and (1.4.21), 


= V ’ “xy = V 

•Scz = -(Ky/Kj.)(8z/8x) Mgy + (3z/8y) 
Similarly, 

®Vy “ 

Myg « (8z/5x) l^y - (VV 


( 1 . 4 . 22 ) 


(1.4.25) 



19 


1,4.3. Stress Resultants for Shallow Shell: 

Considerable simplification of equations (1,4.15) to (1,4.25) 
can be made by assiMing the shell to be shallow |6j. This enables 
us to neglect terms (Qz/^x)^ and the like. Hence, 


and, 


a n d , 


and, 


T 

sx 

5= 

+ (.Wsx) 

T 

sy 


(a^/ey) 

T 

SK 

=S 

-(az/ox) T^^ - (az/ey) 

iV 

ty 



T 

tx 

= 

(dz/ax) Ty^ + Ty^ 



-(az/ay) Tyy -(az/ax) T^y + Ty^ 

^x 

= 

1 

sx 

\y 

xs 

sy 


93 

-(8Z/0X) + (Sz/Sy) Mgy 


« 


V 

ss 

®^x 


=5 

-(8z/6y) M^y + (8z/ax) 


(1.4.24) 


(1.4.25) 


( 1 . 4 , 26 ) 


(1-4.27) 


The equations (1.4.24) to (1.4.27) the forces and 

moments on the shell element in terms of those quantities on 
the plane element. This is done so, because analysis for a 
plane element is much more simpler than that for a curved ele- 
ment. To get the final restate for the shell, the plane, components 
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need these transformations which involve some algebraic manipu- 
lations only. 

Before proceeding further, it will be advantageous if 
the equations to the surface is written out and the necessary 
modifications made. A hyperbolic paraboloid surface can be 
expressed as 


z = (cxy/ab) (1.4.28) 

where, x and y axes represent the characteristic of the hyper- 
bola. Assuming the surface ooasiderably simplieies the procedure 
since in this case we have = «>. Consequently, takiwgcare 

of the cross-symmetry of the shear stresses we have 

= T. = S 
sy tx 

and (1.4.29) 

“sy = “tx = 

1.4.4 Equilibrium Equations: 

The equilibrium e quat ions are written for the plane 
element. 

(6T^/3x) + (eTjj-y/Sy) + ^ 

(dTy^/ax) + ( 93 ?yy/^y) + qy = 0 (1.4.30) 

(dT^yax) + (eTy^Sy) + q^ = 0 

where , 


( 1 . 4 . 31 ) 
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and, 

+ (eM^/3y) - T 3 .Z “ ° 

(1.4.32) 

(6M^/6x) + (aMyy/8x) - = 0 

1.5 IPiPLICifDIOH VABIATIONIL PELWCIPLE 20 HIPAP OHELIS; 

1.5.1 Total Potential: 

To derive appropriate fo2?ce~displacement relations the 
principle of minimum complementary energy with the notion of 
Lagrangian multiplier is used 110|. The functional 

% = (1/2E) / W . dv - /P.U ds (1.5.1) 


is variated subjected to the equilibrium relations to yield the 
appropriate ; .force— displacement relations, where. 


3= <5^ 4- 0^ 2 Ci 




( 1 . 5 . 2 ) 


To obtain an analogue expression for the shell element in terms 
of the pseudo-stress resultants 7,-: consider a rectangular element 
(Pig. 1.8) on the middle surface of the shell. Let 


^12 * ^sy 

M.J, Mg, 1^2 


^21 “ ^tx 


S, 1^3 = T 


se; 


and N 


25 


T 


te: 


and 


^tx ® ^ actual stress resultants 

acting over the faces of the element. 


“21 “ ^sy 


+h/2 I-, 

/ W dv = / / / ((— 

-h/2 


^^)2 + 


h 
12M, 


« )' 


-2 ^ (“^ + — ^ ^ ) Cx — 3^ ^ ^ 


^12 ^21.2 ,2 ^^13 M 

+ (1. )4) lfl+{||)2+(|f)^ ax dy d? (1.5.S) 

4h^ ^ ) 



s 



FiG 1.8 RECTANGULAR & PSEHDO STRESS RESULTANTS ON THE 
ELEMENT OF SHELL SURFACE 
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Equation (1.5.3) can be rewritten, after carrying out the 
integration through the thickness;, as 

/ W dv = 1 // [(N^+Ng)^ + Q '3xdy 

+ ^ // + 2(1 +:^)(m 22-M^M2)J Q dxdy 

+ // + ^23) Q dxdy (1.5.4) 


Where the first double integral on the right hand side denotes 
membrsne energy, the second denotes bending energy and the 
third part is the shear contribution to the total coiapiementary 
energy of the shell element. It is this part which is of direct 
attention to us, about which more is discussed elsewhere in this 
work. 

N. ,Hp etc. are calculated in terms of the pseudo-stress 

I £L 

resultants |7{. 'From Fig, (1.4), we have j?! 


H 

N 

N 

N 


^ Sin w 
2 Sin CO 

15 

25 


= T + S cos CO 
sx 

® ^ COS CO 

ty 



(1.5.5) 


1.5.2; lotal Potential of Shallow Shells: 

Shallow shell approximation on the line of Vlasov enables 
to simplify the energy expressions considerably. A further 
simplification is necessary considering the enormity of the algeb 
raic manipulations demanded for the solution, This constitutes 
in assuming that magnitude of shear stresses are small such that 
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when it is miiltiplied by quantities like (6z/8x), (9z/6y) etc, 

thtt 

iltte resulting product is of higher order. This apparently 
’innocent looking assumption seems to he valid only in cases 
when the shell is ’shallov/’ as well as ’moderately thick’. Though 
no quantitative statement can he made at this stage about the 
natural questions ’how thick’ and ’how shallow’; hut following 
Re issuer jS] it can he said that as long as 9z/6x or 6z/9x does 
not exceed 1/5 and h < 0.05 times the smaller span these assump- 
tions are good. 

However, with this in view equations (1.3.24) and (1.5.25) 


can he 

rewritten 

as 








T = T 

SX XX 

> 

f = T 

sy xy 

and 


dz 

“ - IS 

T 

XX 

- T +T 

oy xy xz 

and 

m 3 - m 

ty ^yy 

f 

ft 

N 

and 


5z 

T 

yx 

— 3} 4"T 

^ yy yz 

and 

^sy ” ^tx 

SB 

T as m 

yx xy 

BS 

S 



(1.5.6) 


with equations (1.5.6) and (1.5.5), keeping in view of the assump- 
tions with reference to equation (1.4.6), 


(Nl + Hg) ■" ^xx ^yy 

Ig = T^^ . Tyy 

N = S 
“l2 

= -(8z/ax) T^ - (92/9y) S + T^jg 

^23 “ -(9z/Sx) S - (0z/Sy) T^y. + Ty^ 
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and. 

Ml Ha = 15^ + Myy 

*12 = “S; (1.5.8) 

*1 • “2 = “« • V 

How, the complementary energy expression (1,4.4) can be written 
as: 

/ W dv = // [(Txx-HDyy)^+2(1+*>)(s2-a}^^ "^yy)]"^^ 

■" 2(l+*>)(Hf-Knc “yy^J ^ '^y 

J'J' [(’^L +Oz/8y)S) 

- aiyg. ((8z/ex) s + (ez/ay)Tyy) ax dy (1.5.9) 

The equation; (1.5.1) is varioted with equation (4.5.9) 
in mind and taking care of the equilibrium equations (1.4.50) and 
(1.4.32) after being multiplied by A’s. Here one point should 
be noted that equation (1,4.32) is not an independent set of 
equations? it can be reduced with the help of (1,4.30c) and hence 
there are three independent equations out of the five as is appa- 
rent £rom (1.4.30) and (1.4.32) and so three values of A, namely 
A.J, Ag SQd are required. The Lagrangian multipliers can 
be distinguished as u, v, w, respectively the three displacement 
functions jlO| for the element. 

Collecting like terms and thin equated to zero, we have 
the force-displacement relations as: 
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(T - T ) - .12(1+^) |z 
^ XX yy-' 5 


5 ^xz 


m, 

Bh ^ 


yy XX 


) , 1 2 . yti>) = Eh|| (1.5.10) 


6 / Sz 


+i2T' 'Is .-,-^-Q . r4iSi + 

5 ^ ^ xz ^ yz"^ 2(1+'^) ^ ox' 

and, (M^- l^y) - 8f^/6x = - ^ d^w/dx^ 


Eh /cm . ^v^ 


(Myy- 1^^) - Jo (1+^) efg/ey^ - ^ aV^y'" (1.5.11) 


Eh' 


v2 6f^ af^ ~3 o 

To ^“5y"*' ■' ^ 


where, 


^1 “ ^xz “■ ^ ^xx " Jj ^ ~ ^ 


dz 


8% 


(1.5.12) 


fg (X, y) = Ty^ - (If) ® - (|f’ ''yy ' 

Eri'ae.tions (1.5.10) and (1.5.11) reduces to that obtained by 

bv Setting 

Reissner for a plate l9{ ^(az/ex), (ez/dy) etc. equal to zero, and 
are obtaim^dL 
boimdary conditions ^as: 

along X = constant 

either, 1 = 0 


XX 

S « 0 


xz 


or u - Ti = 0 

% ^ 

V - T = 0 

w - w = 0 




= 0 


= 0 


a C.M,,vN (1.5.15) 

•. fiC-- =®x 

'♦' ■ ^51^ ^2^^’ “ ®7 



27 


where, 

u, V and w are the translation of the support, 9 and 
0y are the slopes in x and y directions respectively along the 
edge X = constant. !Phe boundary conditions along y = constant 
are similar. 



CHAPOER II 


SOLUTION! 3Y RIT2 TSIROS 


2.1 Ii::.-:CiUCIION: 

It can be shown jlO| that the functions that minimise 
the complementaxy energy as given by equation (1.5.9) also 
determines the solutions of the appropriate Euler's equations, 
in this case the Euler's equations being the compatibility 
equations. Ritz method essentially consists a in finding such 
a function in the form of an infinite series. The advantage 
of this form of the function lies in the fact that the series 

i 

can be terminated after the desired degree of accuracy is 
achieved. The convergence of this sequence to the exact solu- 
tion is usimlly a guaranteed, phenomenon, may be at a cost. 

The method can be described, without going into mathema- 
tical polemics, by considering first the functional to be mini- 
mised as 

I(u) = ff f(x, y, u, u , u , ...) dx dy (2.1.1) 

2 X y 

where 2 spans the whole region under considerations. 

Let the admissible solution u*(x, y) satisfies the 
boundary conditions u = 0(s), on the boundary » and makes I(u) 
a minimum. Row construct. , a sequence 

\ (x, y) - (x, y, a.,, Sg, ..., a^) 


( 2 . 1 . 2 ) 
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with parameters Then construct l(u^) with (2.1.2) and deter- 

mine the parameters a. such that I(u ) is a minimum, the ■Uv,Cx,y) 

* n=0 “ 

forms a solution to the ahove problem. The seq^uence (2.1.2) is 
normally so chosen as to satisfy the prescribed boundaiy condi- 
tions. 

The apparent docility of the method presents a too simple 
outlook, but in practice, the convergence criterion poses a stiff 
obstacle for mnumerical procedure. This is due to the fact 
that the tacit assumption about the completeness |1lj of the 
sequence (2.1.2) is not always possible to achieve in practice; 
However, the loss of accuracy is not alarming compared to that 
already incurred in arriving at the energy expression, and thus 
lends credence for its use here. 

2.2 POTZSTIil. IIERGT IN TERSSS OP DEPORMATIOH: 

The equations (1.5.10) and (1.5.11) form a coupled set 
and hence their solution, if at all possible, involves excessive 
labor and mathematical sophistry. A practical way seems to be 
the minimisation of energy, expressed in terms of displacements. 
Here again, the two membrane displacements, u and v, can not easily 
be separated out finm the displacement in the z-direction (w). 

This difficulty can be obviated if we assume that u and v are 
very small in comparison to w* which is in fact the case for 
shallow shells. With this simplification in mind and the ordi- 
nary plate expressions for moments, we can write equations 
(1.5.9) in three parts as: 


* This is discussed elsewhere in the work. 
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Bending Energy; 


L ex*^ 6y^ 


5^w \2 


^ 2 ( 1 -.) ( 5 ^) 


8^w 8^w ^ 

a 2 * a 2 ^ 

6x 8y^ J 


(ii) oBear Energy: 


u = ff J 

® 120(1- 8x^ 8x 8y^ 


Eh^Kw I dz ( 8:^ e-^w 
5(1-2>^) I 6x^ 6x 0y^ 


/6^w . 8\ n2 . ;r 8^w , 8^w ^2 

+ -^) + ^ — 5 + — » ; 

8x^ 8x By'^ dx'^dy 6y^ 


(iii) stretching Energy (due to w) 
\ “ 11+^ ^ 


)4- || ( + % )| dxd3 

8x^8y 6y^ J 


(2.1.3) 


where. 


o/ah 


(2.1.4) 


(2.1.5) 


Now the functional to he minimised becomes 

X 1 rrF-nL^^w . 5 ^wn2 0/1 ( S^Wn 2 ^fw 

" “ L I I? ^ ^ ■ 8=^2 • 8y2 

8:.5y2^ ^ ey^' f 


,81 f£" + -S^SLl + |2 (— 1 ^ + ^) 

^ dx^ 8x9y^ 8x^9y dy^ 

^ - 2w (Dy^ w - q^) dx dy (2 


( 2 . 1 . 6 ) 
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where, 

— 4 _ „ 6 ^ , 8 ^ 

V “ — I 2 — 5 — s + — T 

0x^ 5x‘^8y^ dy^ 


and, 

B = 


= 



Eh^ 


12(1- 


Ih^d 

+ 1?) 

60 (1 

-^2)2 

Eh^k 


io(.i-" 

51 

Eh k^ 

T* 


2 ( 1 +:^) 


(2.1.7) 



CENTER III 


3.1 SOLUTICI^S TOR PARIICDLAR CASES: 

The displacement function (w) is assumed, satisfying the 
boundary conditions for the case under considerations. With th&s 
displacement function the expression (2.1.6) is written out 
explicitly. The first derivative of the integral, with respect 
to the undetermined constant ia w, is then equated to zero to 
yield the value of the constant for which I is minimum. 


3.1.1 All edges are simply supported: 

The displacement function assumed in this case is the 
usual double sine series: 


00 oo 

m*1 n=1 


Tutting; this value of w in equation (2.1.6) and equating 
to zero, the value of A^ is obtained as 


(3-1.1) 

81 




32 (qo/D)bi__ 
( 2m+1 ) ( 2n.+1 )it 

2_2 




( 2m+1 )^+a^( 2n+1 )^j ^ 


( ((2m+l)^+a2(2n+l)2) -1) 


^ ((2m+l)^ + a^(2n+1)^)) 


( 3 . 1 . 2 ; 


where, 


a « a/b, p = cA» 




(5.1.3) 



a f= Larger Span 

b = Sinaller Span 

c = Mse of tbe shell 

h == Ihickness of the shell 

= Uniformly distributed load 

5.1*2 Two Opposite Edges are Simply Supported and Other Two 
Edges are Pixed: 

The assumed displacement function is 


oo oo 


w =» E E (1-cos Sin 

m=0 n=1 


nTix 

a 


(3.1.4) 


and is oTsrtained as 


2(3. 


b Vi4 




b^ / 


poop 
4a +• 


+ 5aVn^ + 5aVn^ + a^n^) 
a^p^n:^ ((3m^toc^n^) (1+2(-l)^) ) 


5 


(5.1.5) 


- 18 (1- ■»>) ap 

5.1.3; All Edges are Eixed: 

The displacement function is assmed to be 

, = 2 E (1-oosJf2) (5.1.6) 

m=0 n=0 

Proceeding as before, is obtained as; 



V = b4y|£3/(2m+l)^ (1+ ^ + 5aV(2n+1)^ 

4 * 3 

(1+ ( 232+1 ;ii 5 ■*■ 5'n~) ((2nH-l)^+a^(2n+1)^ + 

( 211+1 )^(2m+1)^ ((2m+1)^+a^(2ii+1 )^) ) 

-6ccP^7i:^ (2g (5(2m+1)^ + a^(2n+l)^) + 2((2m+1)^ + 
3a^(2n+1)2)) + 54 (l-.«*>) (3.1.7) 

3.2 ACCUmCT 01 THE SOLUTIONS; 

The accuxacy of the equations (1.5.10) to (1,5.12) were 
discussed before. 

The numerical values obtained fi'om (3.1.1) are resubstituted 
in (1.5.11) and the residue (the diffei^ence between the two sides 
of the equality sign) was calculated as a percentage of the maxi- 
mum moment and was found to be insignificantly small. 

Eor p =» 0 (i.e. the shell reduced to a plate) the moments 
and deflections were checlred with Timoshenko jl2| and G-oldbergl 15| 
as is shown in the accompanied graphs of moment coefficients in 
the central zone (zone I) and are found to be in well ag^ement. 

For a shell with a = 1.5, p = 0.25 and = 0.025 moments 
and deflections were calculated along the centre of the span in 
the shorter direction. These were then compared with that obtained 
by Bouma’s j13j approach and are shown in Fig. 3.1. 

For a shell with » a = 1.0> p = 0.3 and 'i = 0,02, 
w » 3000 Kg, Dayaratnam et al j44i obtains 




2 mm. 


W — " 

max 

Iccoirdiiig to the pi^esent inTestigation 


^max ~ ^ 




^ax 


* 4» 15 Kg. cm. 


1.3.3 Discussion of the Design Graphs: 


Three aspect ratios have been defined as 

a - a/b , the Span ratio 

P = c/b , the rise ratio 

T? = h/b , the thickness 3atio. 


Moments and deflection coefficients have been plotted for 
various values of a, p, . These deflection coefficients refer 
to the deflection of the middle point of the shell where the 
maximum deflection occurs. These coefficients are non-dimensi- 
onalised parameters and the exact relations for moments and def- 
lections are, 



6 (q/E) 1 

(3.3.1) 


1 1 

(3.3.2) 

ly " 

^ ly 

(3.3.3) 


The subscripts in equations (3.3.2) and (3.3.3) refers to the 
zone nxjmber and direction of the vector respectively. To facili- 
tate the design and placing of reinforcements, the shell is 



divided into nine aones (Pig. 1.4) of which four are typical. 

Por each value of span ratio, the curves are the variation 
in moment coefficients for different values of thickness ratios 
i '^ ) and for a particular zone and direction, the corresponding 
rise ratios (p ) are also specified therein, 

The apparent strangeness in presenting only the bending 
moments without a mention of the membrane forces can be explained 
on two grounds: 

1* The effect of membrane forces in a thick shell is negli- 
gible in the sense that the minimum allowable reinforcements (as 
is prescribed by IS codey far exceeds that obtained from 

an actual membrane stress calculation. 

2, Prom an investigation ll7| it is inferred that for p 

ratio lying within 0.025 to 0.001 and p varying between 0.1 and 
0.25 the moments predominates to a high degree. The cases inves- 
tigated in this work fall in this category and thus substantiates 
the prescription of a design method based on bending moments. 

3.4 IKBuiaN PROCEDUEE: 

The design procedure can best be illustrated by an sample. 
The superiority of the shell fora dat ion can be seen from the com- 
parison of the two designs presented below, the first one is for 
that of a conventional footing and the second one is the shell 
foundation which is advocated in this work. 

Example: Suppose a footing has to be designed for a 

column carrying a load of 25 tonnes in a 


25 cm X 25 cm. square 



soil of bearing capacity of 1.0 T / (or 0.1 Kg/cm^). 
Conventional design: (Mg. 5,35) 

Suferln^osed load = 25 tonnes 

Footing load (10 percent) = 2.5 Tonnes 

Total vertical load 27.5 Tonnes 

Area required = 27.5/1.0 = 27.5 

Provide: 5.1“^ x 5.5“ = 30.25 ^ 

Actual soil press-ure = 27.5/30.25 = 0.91 T/M^ <1 T/M^ 

(Safe). 

Max.“ bending moment = (I/24) x 0.91 (2 x 5.5 + 0.25) 

(5.5 - 0.25)^ 

= 11.75 X 10® Kg cm 

Depth required = ifCH. 75x10 )/(25x9.7) = 48.5 cm. 

Depth provided at the base of the column = 50 cm, 

(See Ilg,3.36) 

Intensity of Punching shear = (1/4) x ^ * ^0. 15x0,185^^ ^ 

= 5.66 Kg/cm^ < 10.5 Kg/cm^ 

(Safe) 


Depth provided at the edge = 15 (See p^ig. 3.36) 

Effective depth at the critical section = 35.7 cm. 

Shear force = ?10 x 6.p x 2,1 2 5 „ 5520 Kg 

Shear stress = igs^xSO" ‘ 



Area of steel 


52.5 Sq., cm. 


11.75 X 10 - X 1 

1265 X 6 X 48.5 

Use 16 mm ^ at the rate of 20 cm. c/^ . (See Pig. 3.56) 

Shell Design: 

The area of the shell in plan is assiamed to he the same, 
since this is governed by the allowable bearing pressure, 

Jbr this case, 

Span Batio, 

a = 5. 5/5. 5 = 1.0 

a uniform shell thickness of half that of the ’plate* is assumed 
So thickness provided = 25 cm, 

thickness ratio, = 25/550 - 0,05 

Pigs. 3.8 - 3.10 and Pig, 3.4 gives the moment and deflection 
coefficients for the span ratio of 1,0. 

Though a greater rise will reduce the bending moments 
considerably but to show the effectiveness of the method only 
10 percent rise (i.e, p = 0.1) is assumed. 

Prom the above graphs, and for these values of p and ^ , 
the following values are obtained 

6 * 300ii = 5.23, = 5.52, Hi * 4-.94 

Corresponding moments: 

% = 15.90 X105 Kg cm 

= 16.70 X 10^ Kg cm 
15.00 X 10^ Kg cm 


Mjii= 



Maximum depth required = 8 cm < 20 cm (Safe; 

a. 

But SO smaJ.1 ^depth is not feasible considering the punch— 
ing shear. With the assumed depth of 20 cm, the intensity of 
punching shear is obtained as 11.66 Kg/cm^ ( 10.5 Kg/cm^ and 

hence safe). 

Area of steel required will be very small considering the 
ratio of maximum bending moments in the his cases. Hence, the 
mtoimum permissible steel is to be provided (i.e, 10 mm at the 
rate of 25 cm. ^/e). 

Abave two cases clearly show the superiority of shell 
foundation over ordinary spread foundation, 

3*5 DISCUSSION 01' EE3IILTS: 

The deflection coefficients (&) decreases very rapidly 
as the thickness ratifli ( ) increases from 2.5 percent to 5 per- 

cent (Pig. 5.4 - 5.7). This decreement rate then becomes asympto- 
tic for further increase of thickness ratio. Por 10 percent rise 
of the shell (i.e. p =0,1) the deflection coefficient reduces 
to 1/5 to 1/8 for thickness ratio varying between the above two 
limits; the later values are for higher span ratios (i.e, oC > 1.0), 

With only a 10 percent rise the shell deflection coeffi- 
cient reduces considerably from the plate values (See Pig. 5.6 
and 3.7). Increasing the shell rise also reduces the deflection 

* Punching shear in a shell will be smaller due to the vertical 
component of the inclined membrane forces, but this effect was 
not considered here for simplicity. 



coefficient (6) considerably for tbe thin shell range. Eor 25 
percent shell rise the 5 Talues are almost half of the 6 values 
for 10 percent rise. But the decreements are not so much observed 
for higher thiclcnesses; in fact, for 7.5 percent or more thickness 
ratio the shell rise has no perceptible effect on deflection co- 
efficients, This later observation can be explained from the 
point of view that with increase percentage thickness the curva^ 
ture effect of the shell decreases and hence its ‘shell superiority 
over the plate is decreased reducing it almost to an inclined plate 

That with increased thickness the shell starts ^‘behave - al- 
most as a plate, particularly for small rise to span ratio, is 
also corroborated by the dairies for moment coefficients. For 
example, all the curves show little perceptible difference between 
plate moment values and that for shells near the 10 percent thick- 
ness range, though for thin shell values a wide difference is 
observed, as is expected. 

A look at Fig. 5*8 to 3.34 shows that a 10 percent rise 
reduces the moment coefficients by at least 30 percent from the 
plate values. This reduction normally lies between 30 percent 
to 50 percent. With increase in shell rise the moment coeffi- 
cients reduces considerably; this effect is particularly pro- 
nounced near the central zone (zone I, Fig. 1,4). 

A con^arison of Fig. 3.11 and Fig. 3.15 shows that increas- 
ing the span ratios ( ) decreases the moment in the longer direc- 
tion. This effect is particularly pronounced for a > 1.5 and near 
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the central aone (aone I, i*ig. 1.4.), 

One particular case of a shell with a fixed - fixed 
edge condition was ntnaericaliy solved. !I?he shell is of 4M x 6M 
in plwa. fhe bending moment variation along the centre of the 
longer span is shown in Pig. 3.5. It is observed that the edge 
moment varies from 110 percent to 140 percent of the central 
moment for shell thickness ratio f varying from 2.5 percent 
to 10 percent and for a shell rise of 25 percent. 

3.$ COIJCL’JSIOJ; 

Since the assumption of simply supported boundary condi- 
tions can not exactly be realised in practice, this inherent 
limitation of the graphs restricts its use for design office 
practice. Nevertheless, this work points out several facts 
conclusively; 

(i) The reduction in volume of concrete and steel is consi- 
derable, as the design example also illustrates. 

(ii) Even for foundation the shell thickness is restricted 
more by practical considerations than by bending stresses, 
as is usual with thin shell design practice. 

(iii) The area of steel is governed by the allowable requirement 
of the code of practice rather than the stresses occuring within, 

(iv) Nothing has been said :.about secondary effects such as 
settlement, hydrostatic lift etc. which are vitally important in 
practice. So befo2?e its adoptation for design purpose a much 
more thorough and detailed investigation is suggested. 



fliia worfc point© ont to the Justification for adopting 
thin shell design method for foundation shells also, as is usually 
done* Only difference lies in the considerations of bending 
moments near the comers which may at times be high. Ihis once 
again illuminates . the path to be followed in future work. 



HABIiE 1 


Bellection and Hoisent Coefficients for Simply supported Hyper shell. 


p 

i 

6 


/"ll 


/•it 

§.1 

0.025 

1640 

4.37 

4.35 

3.30 

4.35 


0.050 

3.00 

5.23 

5.52 

4.94 

5.52 


0.075 

108 

5.63 

5.96 

5.46 

5.96 


0.10 

48 

6.01 

6.23 

5.68 

6.25 

0.15 

0.025 

1310 

4.06 

3.88 

2.69 

3.88 


0.050 

280 

5.06 

5.26 

4.55 

5.26 


0.075 

100 

5.63 

5.80 

5.24 

5.80 


0.10 

48 

5.95 

6.14 

5.54 

6 , 14 

0.20 

0.025 

1150 

3.83 

3.55 

2.27 

3.55 


0.050 

280 

4.91 

5.03 

4.23 

5.03 


0.075 

102 

5.43 

5.67 

5.03 

5.67 


0.10 

44 

*5.89 

6.04 

5.40 

6.04 

0.25 

0.025 

960 

3.64 

5.29- 

1.96 

3.29 


0.050 

260 

4.77 

4.83 

3.94 

4.83 


0.075 

98 

5.35 

5.54 

4.84 

5.54 


0.10 

43 

5.82 

5.96 

5.27 

5.96 


—2 

ill yM -values should he multiplied by 10 . 


lote : 


’lABIE 2 

deflection end toent Ooeffioiente for elnply Supported Hyper ebell 


p 


6 

Ax 

llz 

^ IVx 

^ Ey 

u 

' llj 

'^Illy 

iVy 

0.0 

0.025 

1568 

2.82 

3.66 

3.10 

2.09 

4.45 

3.82 

3.97 

3.96 


0.050 

172 

2.82 

3.67 

3.13 

2.10 

4.f6 

3.83 

4.02 

3.99 


0.075 

51 

2.83 

3.73 

3.20 

2.11 

4.48 

3.85 

4.12 

4.04 


0,10 

22 

2.83 

3.78 

3.30 

2. 12 

4.49 

3.86 

4.26 

4.12 

0. 10 

0.025 

191 

1.61 

2.74 

2.50 

1.32 

2.66 

2.59 

3.17 

2.80 


0.050 

141 

2.38 

3.36 

2.93 

1.82 

3.80 

3.39 

3.74 

3.57 


0.075 

46 

2.61 

3.56 

3.10 

1.97 

4.15 

3.63. 

3.97 

3.83 


0.10 

20 

2.70 

3.68 

3.24 

2.04 

4.30 

3.73 

4.17 

3.99 

0.15 

0.025 

661 

1.30 

2.48 

2.33 

1.15 

2.21 

2.27 

2.94 

2.50 


0.050 

131 

2.20 

3.23 

2.84 

1.71 

5.55 

3.22 

3.62 

3.40 


0.075 

44 

2.51 

3.49 

3.05 

1 .91 

4.00 

3.53 

5.91 

3.74 


0.10 

19.7 

2.64 

3.64 

5.21 

2.00 

4.21 

3.67 

4.13 

3.93 

0.20 

0.023 

570 

1.09 

2.28 

2.19 

0.98 

1.90 

2.03 

2.77 

2.28 


0.050 

127 

2.05 

3.11 

2.77 

1.61 

5.32 

3.06 

5.52 

3-25 


0.075 

43 

2.42 

3.42 

3.00 

1.85 

3.87 

3.44 

5.85 

3.65 


0.10 

19.3 

2.58 

3.59 

5.18 

1.96 

4.12 

3.60 

4.09 

3.87 

0.25 

0.025 

484 

0.92 

2.11 

2.08 

0.87 

1 .66 

1.83 

2.63 

2.11 


0.05 

116 

1.92 

3.01 

2.70 

1.53 

3.12 

2.92 

3.43 

3-13 


0.075 

42 

2.53 

5.36 

2.97 

1.80 

3.74 

3.35 

2.79 

3.57 


0.10 

19 

2.52 

3.55 

3.15 

1.92 

4.03 

3.55 

4.05 

3.81 


Hote : All values should be multiplied by 10”^ 



giBIiE 5 

Deflection and Moment coefficients for simply supported hyper shell. 


ij 

P 


6 



^ riix 

'lYx 

'^ly 

^Ily 

^iiiy 

^lYy 

0.0 

0.025 

568 

1.48 

2.36 

1.95 

1.09 

3.17 

3.03 

2.92 

2.65 


0.050 

71 

1.48 

2.37 

1.95 

1.10 

3.18 

3.04 

2.95 

2.66 


0.075 

21 

1.48 

2.39 

1.98 

1,10 

3.19 

3.05 

3.01 

2.69 


0.10 

9 

1.49 

2.42 

2.03 

1.11 

3.20 

3.06 

3.09 

2.73 

0.10 

0.025 

368 

0.90 

1.86 

1.61 

0.75 

2.08 

2.23 

2.41 

1.95 


0.050 

62 

1.29 

2.21 

1.84 

0.97 

2.81 

2.77 

2.78 

2.43 


0*075 

20 

1.39 

2.31 

1.93 

1.04 

3.00 

2.92 

2.93 

2.58 


0.10 

8.6 

1.43 

2.37 

2.00 

1.07 

3.09 

2.99 

3.04 

2.67 

0.15 

0.025 

315 

0.74 

1.70 

1.51 

0.63 

1.78 

1.99 

2.25 

1.75 


0.050 

60 

1.20 

2.14 

1.80 

0.92 

2.65 

2.66 

2.71 

2.33 


0.075 

19 

1.34 

2.28 

1.91 

1.01 

2.93 

2.86 

2.89 

2.52 


0.10 

8.5 

1.40 

2.35 

1.98 

1.05 

3.04 

2.95 

3.02 

2.63 

0.20 

9.025 

277 

0.62 

1.57 

1.42 

0.55 

1.55 

1.81 

2.12 

1.60 


0.050 

56 

1.13 

2.08 

1.76 

0.87 

2.51 

2.56 

2.65 

2.24 


0.075 

19 

1.30 

2.24 

1.88 

0.99 

2.85 

2.81 

2.85 

2.48 


0.10 

8.4 

1.38 

2.33 

1.97 

1.04 

2.99 

2.91 

3.00 

2.60 

0.25 

0.025 

246 

0.53 

1.47 

1.35 

0.49 

1.37 

1.66 

2.02 

1.48 


0,050 

53 

1.06 

2.02 

1.72 

0.83 

2.39 

2.47 

2.59 

2.16 


0.075 

18 

1.26 

2.21 

1.86 

0.96 

2.77 

2.75 

2.82 

2.43 


0.10 

8.5 

1.35 

2.30 

1.96 

1.02 

2.93 

2.88 

2.98 

2.57 


Mote: ill -values should he multiplied hy 10' 
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Deflection and Moment Coefficients for Simply Supported Hypar Shell. 

a * 2.5 


p 


5 

^ Ix 

/u 

• IIx 

/^IIX 

^Yx 

/ly 


Clly 

h 

lYy 

0,0 

0.025 

269 

0.86 

1.62 

1,30 

0.64 

2.28 

2.42 

2.22 

1.84 


0,030 

34 

0.86 

1.63 

1.31 

0.64 

2.29 

2.42 

2.24 

1.85 


0.075 

10 

0.86 

1.64 

1.33 

0.64 

2.29 

2.43 

2.28 

1.87 


0.10 

43 

0.86 

1.65 

1.35 

0.65 

2.30 

2.44 

2.34 

1.90 

0.10 

0.025 

192 

0.56 

1.33 

1.11 

0.45 

1.62 

1.90 

1.89 

1.43 


0.050 

31 

0.76 

1.54 

1.25 

0.58 

2.08 

2.26 

2.14 

1.72 


0.075 

9.7 

0.81 

1.59 

1.30 

0.62 

2.19 

2.35 

2.23 

1.81 


0.10 

4.2 

0.83 

1.63 

1.34 

0.65 

2.24 

2.39 

2.31 

1.86 

0.15 

0.025 

169 

0.47 

1,23 

1.05 

0.40 

1.42 

1.72 

1.78 

1.30 


0.050 

29.7 

0.72 

1 .50 

1.23 

0.56 

1.99 

2.19 

2. 10 

1.66 


0.075 

9.4 

0.79 

1.57 

1.29 

0.60 

2.15 

2.32 

2.21 

1.78 


0.10 

4.1 

0.82 

1.61 

1.33 

0.65 

2.22 

2.37 

2.29 

1.84 

0.20 

0.025 

146 

0.40 

1.14 

0.99 

0.35 

1.26 

1.58 

1.69 

1.19 


0.050 

27.8 

0.69 

1.46 

1.20 

0.55 

1.90 

2.12 

2.05 

1.61 


0.075 

9.1 

0.77 

1.56 

1.28 

0.59 

2.10 

2.28 

2.19 

1.75 


0.10 

4.4 

0.81 

1.60 

1-32 

0.62 

2. 19 

2.35 

2.28 

1.82 

0.25 

0.025 

131 

0.34 

1.07 

0.95 

0.51 

1.13 

1.46 

1.61 

1.11 


0.050 

27.0 

0.65 

1.43 

1.18 

0.51 

1.82 

2.06 

2.01 

1.56 


9.075 

9.1 

0.75 

1.54 

1.26 

0.58 

2.06 

2.25 

2.16 

1.72 


0,10 

4.0 

0.796 

1.59 

1.32 

0.61 

2.16 

2.33 

2.27 

1.81 


— — — ^2 

Sote: 111 ^ values should be multiplied by 10 . 



X, y, « 

8,t 


:T02:Arici\r 


i, 3, k 

m 


r 

m 

h 


h 




c 

a 

b 

a 

'i 

E 

W 


= Crthogonal cartesian coordinate system. 

« Surface coordinate system the x and y direction 
respectively. 

= Unit vectors in the cartesian system. 

“ Unit vectors in s, t and normal to the surface 
respectively. 

= Radius vector of a point from a chosen origin, 

= Included angle between s and t direction. 

= Radius of curvature along s and t direction 
3?espectively. 

= Thickness of the shell. 

= Membrane force in i-plane j - direction, 

= Moment vector in i-plane ^-direction following a 
right handed cork-screw rule. 

= Force in i-plane 
= Moment vector in i-plane 
= Rise of the shell. 

= Span in x direction (longer) 

= iSpan in y direction (Smaller) 

= a/b = Span ratio 

= c/b = Rise ratio 

= hA> = Thickness ratio 

= Modulus of elasticity of the material. 

= Complementary energy per unit volume. 



o*^,cr 2 

z 

i3 


Normal stresses in 1 and 2 direction. 
Shear stresses in i-plane, ;5 - direction. 


11, V, w 

'j> 

0 

D 


* Displacements in x, y and z - direction^ 
» Poisson *s ratio. 

“ Strain energy per unit volume. 

* ilerural rigidity of a plate element. 

== External load vector 


iny other notation used at any phase is defined 
before it is used. 
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111 til® Vais Tolw® tli« relation# haT® 

ti«#a asfuiBed to write the eaerisr equation in teiw# of lii^iaee- 
meiite. Thl# can easily he lezl'wel by writing down the malt 
vectors in g* t and direetion in the defowiad gtato and then 
calculating the strains. 

Ihe surface is given by 

z » zUt y) (At) 

Before defoneatloni a point in the shell can be denoted 
veotorially as: 

r » xl *► 3 ^ + e{x, y) 1 (A2) 

Therefore, the unit vectors in s, t and ^ direction are obtained, 
in undeformed state, as (see equation 1,4. 5). 

I + «_ 1 

•• ■ 777^^ 

^ J + B 1 c 

■ 77TT5^ 

m 

1 ~ z^T ~ 

( 1 +*| ) 1 +sp 

SB- 8J- 

X y 

«■ 0 (A4) 

xt 0 


and 


V . • “y / 


e^ » e X e. 

St 


cos oi 


's • ®t 


cos 0 « 

0 8 





wimmf 

w * angle between s and t dineetlon befert defor^^atlon 
0 • %/Z • angle between s ant t direetioa before 

defomstion 

0 «* Js/2 »■ Mgle between t and % dlreetion before 

defonnation 


After defoa®atlon tbe sane ouantitiee beoene 

r* » r -ef « [x^fuCx.y)] X -f [y+'r(x,y)J J 4- [B(x,f)4 w(x,y^I 

(A5) 

Hlberefore, tbe tmit vectors, In the defoxraed state, ares 

% * / I I 

(t -«> ||)I ^ If 7 ♦ t |g) X 

•; - ifV-iif’i 


1 ^ I t (l+ J 4 (My 4 I®) 1 



Ci0) 


1/2 


* The dleplaceaeate are functions of x, y and s. Mt they 
can easily be reduced to functions of at and y for tbe oaee 
under consideration with the help of equation (A1), 



m 




jS* / I 
is is 



She equatloa (A6) gets nmch more stolified after tke aoiwliaoar 
tens are eliialaated and usual shallow shell approximations aue 
mMe* With these in mind (A6) reduces to 




?4 


» gg ^ 


^ If + *y ^ 


Co8 9’ » 2^ (*1^ 4- 4. « Jsi 


(*if + + By (U2 |S) 

(1 ♦ 2 |2) + *y (|B + ^) 


Ulo) 


How* for a shallow shell* 'the stzialn coniponent s are ehtaiaed 


(equating e to zei:o)i 

sum 


€ mS^ 

yy ©y 


0u . 

^5 ♦ *: 


^ * *y ^ 


(Atl) 


Bjt ^ly ^ ®y If) 


(tig) 


Vb “ ♦ »y Clf ■»‘^) 


JPhese are the stsraiz^dlsplaeeraeiit relation in this 
work to express energy in terra© of dtsplaowents,. 



